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1. Introduction 

The doubling problem of fermions has been one of the notorious obstacles to QCD simulations. 
By now several fermion constructions to bypass it have been developed, although all of them have 
their individual shortcomings. There is also an alternative possibility to settle the problem. More 
than 20 years ago, Karsten and Wilczek proposed the minimally doubled action [|I|], where only 
two species emerge. Recently one of the authors of the present paper [Q] has also proposed a 
two-parameter class of fermion actions. These minimally doubled fermions all possess one exact 
chiral symmetry and exact locality. As such they should be faster for simulation, at least for two- 
flavor QCD, than other chirally symmetric lattice fermions. In this paper we propose a systematic 
method to reduce the number of doublers, which we term twisted-ordering method [|p. In this 
way we obtain new classes of minimally-doubled fermions. We classify all the known minimally 
doubled actions into two types. By this classification we can derive several further minimally 
doubled actions deductively. We also study higher-dimensional extension of them and find the 
parameter range for minimal-doubling gets narrower with the dimension in most cases 

2. Twisted-ordering Method 

In this section we propose a systematic way of controlling the number of species of lattice 
fermions within the requirement of Nielsen-Ninomiya's no-go theorem [Qj. We will first discuss 
the 2-dimensional case to give an intuitive understanding of this mechanism. Then we will go on to 
the 4-dimensional case and show we can construct a new minimally doubled action by this method. 

Let us begin with the following simple 2-dimensional Dirac operator in momentum space with 
0(a) Wilson-like terms. 

D(p) = (sinpi +cospi - 

+(sin/7 2 + cos p 2 - l)//2, (2.1) 

where p\, pi and J\, 72 stand for 2-dimensional momenta and Gamma matrices respectively. Here 
the deviation from the usual Wilson action is that the 0(a) terms are accompanied by Gamma ma- 
trices. There are four zeros of this Dirac operator at (pi,p2) = (0,0), (0, 71 / 2), (n/2,0), (tz/2,%/2). 
This means that the number of the species is the same as for naive fermions. 

Next we "twist" the order of the 0(a) terms, or equivalently, permute cos p\ and cos P2 to give 

D(p) = (sinpi + cos/>2 - 

+(sin/7 2 + cos/7i - l)/y 2 . (2.2) 

Here only two of the zeros ((0,0) and (n/2,7i/2)) remain and the other two are eliminated. Thus 
the number of species becomes two, the minimal number required by the no-go theorem. As seen 
from this example, twisting of the order of 0(a) terms reduces the number of species. We call 
this method "twisted-ordering" in the rest of this paper. We depict the appearance of zeros in the 
Brillouin zone for the above two cases in Fig. |] and ||[ 

Now we consider excitations at the zeros. The Dirac operator is expanded about the zeros as 

D {l) (q) = iYiqi + ifiqi + 0(q 2 ), (2.3) 
D {2) (q) = -m<?2 - iftqi + 0(q 2 ). (2.4) 
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Figure 1: Red-dotted and blue-solid lines stand for 
zeros of the coefficients of j\ and 72 respectively. 
Black points stand for the four zeros of the Dirac op- 
erator. 



Figure 2: Red-dotted and blue-solid curves stand for 
zeros of the coefficients of j\ and 72. Black points 
stand for two zeros (0,0) {n/2, n/2) of the Dirac op- 
erator. 



Here we expand with respect to defined as p„ = p^ + and we denote the two expansions 
as for the zero (0,0) and for the other at (n/2, n/2). Momentum bases at the zero 
(0,0) are given by b[ = (1,0) and b 2 = (0, 1) while those at the other (n/2, n/2) are given by 

(2) (2) 

bj = (0,-1) and b 2 = (—1,0). This means excitations from the two zeros describe physical 
fermions on the orthogonal lattice. The gamma-5 at p = (0,0) is given by 7 5 =7x72 while the 
gamma-5 at p = (n/2, n/2) is given by yf^ = 7271 = —75 This sign change between the species 
is a typical relation between two species with minimal doubling since Nielsen-Ninomiya's no-go 
theorem requires fermion pairs to possess chiral charges with opposite signs. 

We also obtain 4d minimally doubled actions with this method. The simplest form of 4d 
twisted-ordering actions is given by 

D(p) = (sinpi + cos P2 - l)z'7i 
+(sinp 2 + cos/73 - 1) i'72 
+(sinp 3 + cos p 4 - l)/73 

+(sin/? 4 + cos/?i - l)/7 4 , (2.5) 

with 0(a) Wilson-like terms whose order is twisted. Here exist only 2 zeros 

(PhP2,p3,p4) = (0,0,0,0), (n/2,n/2,n/2,n/2). (2.6) 

We have confirmed there are no other real zeros of this operator numerically. This is a new type 
of minimally doubled fermion on the orthogonal lattice. In this case excitations from the two 
zeros describe physical fermions on the orthogonal lattice again. And the two associated 75 's have 

v ■ (2) (1) 

opposite signs as 7 5 = 74717273 = —75 . 

Gauging these theories is straightforward. One merely inserts the gauge fields as link oper- 
ators in the hopping terms for the action in position space. Specifically, for the twisted-ordering 
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minimally doubled fermion in position space we have 
1 
2 



(2.7) 

where we define 

3. Minimally-Doubled Fermions 

In this section we discuss two classes of minimally doubled actions obtained from the original 
twisted-ordering action ( |2.5[ ) [Qj. The first one, which we call "dropped twisted-ordering action", is 



constructed as following: We drop one of cos p^ — 1 terms in the Dirac operator (|2.5[), for example, 
drop cos p\ — 1 as 

D(p) = (sinpi + cos/72 - 1) ij\ 
+(sin/? 2 + cos ^ 3 - 1) iy 2 
+(sinp 3 + cosp 4 - 1)1/3 

+ sinp4/y4. (3.1) 

This has two zeros given by 

(pi,P2,P3,P4) = (0,0,0,0), (71,0,0,0). (3.2) 
It is obviously a minimally doubled action. 



Next we consider a second variety where we turn on a parameter in the action ( |2.5[ ) as follow- 
ing, 

D(p) = (sin/?i + cos/?2 — oc) iji 
+(sin/? 2 + cos/73 - a)iy 2 
+(sin P3 + cos />4 — a) ijs 

+(sin/> 4 + cos/>i — a) ify. (3.3) 

Here we replace unity with a positive parameter a in the constant terms of the operator. The 
parameter range for minimal-doubling is given by 

< a < V2. (3.4) 

In this range two zeros of the Dirac operator, which we denote as py and py, are given by 



pP = arcsin I — — I = arcsinA, (3.5) 



a 



V2-a 



a + V2-a 2 



p^ = arcsin — V ^ —J = arcsin \/l — A 2 . (3.6) 

where for convenience we define a parameter A as A = (a — y/2— a 2 ) /2. These two zeros reduce to 



Eq. (2.6) with a = 1. By expanding the Dirac operator we see momentum bases are non-orthogonal 



for cases of a / 1; thus, the associated lattices are non-orthogonal as with the Creutz fermionri 
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4. Symmetries 



In the previous section we have studied two classes of minimally doubled actions in ( |3.1| ) and 



(p.3|), both of which are obtained from the original twisted-ordering action in (|2.5|). In this section 
we will discuss discrete symmetries of these actions and redundant operators generated by loop 
corrections. 

Firstly we note both of the actions possess some common properties with other minimally dou- 
bled fermions: "Gamma-5 hermiticity", "Discrete translation invariance", "flavor-singlet U(l)y" 
and "flavor-nonsinglet U(1)a-" The last one is the exact chiral symmetry preventing additive mass 
renormalization for the neutral pion. Besides them, they have exact locality and gauge invariance 
when gauged by link variables. We expect every sensible minimally doubled fermion to possess 
these basic properties. 

On the other hand, Discrete symmetries associated with permutation of the axes and C, P or 
T invariance depends on a class of the actions. In other words, we can identify minimally dou- 
bled actions by these discrete symmetries. Here we show discrete symmetries which the dropped 
twisted-ordering action possesses: l.CP, 2.T, 3.Z2 associated with two zeros. Here we take the P4 
direction as time. We also write the symmetries which the full twisted-ordering action possesses: 
l.CPT , 2.Z4 associated with axes permutation, 3.Z2 associated with two zeros. Because of lack of 
sufficient discrete symmetries, redundant operators can be generated radiatively in these action as 
in the other classes [Sj. For example, the case of the twisted-ordering action is similar to that of 
Borici-Creutz action which has S4 symmetry [S]. We speculate the number of redundant operators 
generated in the twisted-ordering action is the same as those of the Borici-Creutz one shown in 



5. Classification 

In this section we discuss a classification of the known minimally doubled fermions. So far we 
have seen four variations on minimally doubled fermions: Karsten-Wilczek [jl]], Borici-Creutz 
twisted-ordering and dropped twisted-ordering fermions Here we classify these actions into 
two types: One of them, which includes Karsten-Wilczek and dropped twisted-ordering actions, is 
given by 

where we take i = 1,2,3 as i while j = 2,3,4 as j. The point is that indices i and j are staggered. 
The different actions in this class depend on the choice for the matrix R. For example, consider the 
following R's 
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(5.2) 



For the first R the general form (5.1) reduces to the Karsten-Wilczek action while it reduces to 
the dropped twisted-ordering action for the second one. The action for the third case is a new 
possibility. 
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We see that there are many options associated with possible R's. One common property with 
fermions obtained in this way is that two zeros are separated along a single lattice axis and given 



by (0,0,0,0) and (71,0,0,0). One significance about the general form ( |5.1| ) is that a coefficient of 
at least one Gamma matrix has no (cos — 1) term. It is also notable that, in a coefficient of each 
gamma matrix, a momentum component associated with a sin term differs from the component 
in any cos term. These two points seem to be essential to minimal doubling. 

The second class of actions includes the twisted-ordering fermion with the a parameter and 
the Borici-Creutz actions. A generalized Dirac operator for this type is given by 



D(p) = /£[7 /i sin( j p M + J 8 M ) 
I 1 



Y sm(p^-^)} - iF 



(5.3) 



where y 1 ^ = A^ v y v is another set of gamma matrices where we define A as an orthogonal matrix 
with some conditions: At least one eigenvalue of A should be 1 and all four components of the 
associated eigenvector should have non-zero values. Here and T has a relation with this A as 
r = Yy. sin 2/3^ = T^j sin 2/3^, which means sin 2/3^ is an eigenvector of A as A^ v sin2/3 v = 
sin 2/3^. Thus once A is fixed, /3^ and T are determined up to a overall factor of sin 2/3^. By 
imposing these conditions on A, and T, the action ( |5.3[ ) can be a minimally doubled action. In 
such a case indicates locations of two zeros as p^ = i/3^. Adjusting j8, we can control the 
locations of two zeros. Note that the first and second terms in ( |5.3[ ) are nothing but naive fermion 
actions. We can eliminate some species by combining two naive actions with different zeros in one 
action. 



Now let us show this action includes minimally doubled actions. ( p.3| ) reduces to the Borici 
action and the twisted-ordering action by choosing A's as 





— A^vYv, 
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and 



/0 l\ 
10 
10 

\0 1 0/ 



(5.4) 



where we fix /3^ = 7i/4 and T is given by T = y^. Finally we have classified all the known 
minimally doubled fermions into the two types. Now we can derive a lot of varieties from these 



general forms deductively as we have shown in (5.2) 



6. Higher-dimensional Extensions 

In this section we discuss higher-dimensional extensions of minimally doubled actions. As two 
of the present authors showed in [Q], we can extend minimally doubled fermions to general dimen- 
sions. The main observation in this generalization is that the parameter range for minimal-doubling 
tends to get narrower in higher dimensions. For example, one of J-dimensional generalizations of 
Karsten-Wilczek fermion is given by 

d 

Dkw =£[/y Ai sin/^ + iy d (r- cos p^)}, (6.1) 
l 1 
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where r is a parameter which will be fixed to realize minimal-doubling. Here the parameter range 
for minimal-doubling is given by 

d + V2-3 d + V2-\ 

1 <r < 1 • ( 6 - 2 ) 

d d 

It is obvious that this range gets narrower with the dimension (d) getting larger. This is also the 
case with even-dimensional generalization of Borici-Creutz action, which has two parameters B 
and C. The minimal-doubling range for C is given by 

d-2 



, < C < 1 (6.3) 
d 

where we concentrate only on the case d = 2m (m = 1,2,3...). Here this range gets narrower with 
the dimension. As seen from these two examples, it is more difficult to realize minimal-doubling 
in higher dimensions. It is natural since the number of species increases with the dimension and 
we need to impose a stronger condition to obtain only two species in higher dimensions. On the 
other hand things are a little different in the twisted-ordering fermion. In this case the minimal- 
doubling parameter range for a in Eq. (3.4) holds even in higher dimensions [|3j]. This is because 
the mechanism "Twisted-ordering" itself reduces the number of species in this case while adjusting 
the parameter controls them in the other cases. This fact may imply the twisted-ordering fermion 
is the most fundamental type among all the known minimally-doubled fermions. 



7. Conclusion 

In this paper we propose a method to control the number of species of lattice fermions and 
obtain new classes of minimally-doubled fermions. We also classify all the known minimally dou- 
bled actions into two types. By this classification we can derive several further minimally doubled 
actions deductively. We also study a higher-dimensional extension of them and find a parameter 
range for minimal-doubling gets narrower with the dimension in most cases except for the twisted- 
ordering fermion. With several varieties of minimally doubled fermion actions available, a goal is 
to apply these actions to numerical simulations and study their relative advantages. Such studies [^] 
so far focused only on Karsten-Wilczek and Borici-Creutz fermions. One can now explore other 
varieties such as ones in this paper. 
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